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Abstract. A systematic procedure is proposed for deriving all the gauge symmetries of the gen- 
eral, not necessarily variational, equations of motion. For the variational equations, this procedure 
reduces to the Dirac-Bergmann algorithm for the constrained Hamiltonian systems with certain 
extension: it remains applicable beyond the scope of Dirac's conjecture. Even though no pair- 
ing exists between the constraints and the gauge symmetry generators in general non-variational 
dynamics, certain counterparts still can be identified of the first- and second-class constraints with- 
out appealing to any Poisson structure. It is shown that the general local gauge dynamics can be 
equivalently reformulated in an involutive normal form. The last form of dynamics always admits 
the BRST embedding, which does not require the classical equations to follow from any variational 
principle. 



1. Introduction 

By a local dynamical system we understand the dynamics whose true trajectories are defined 
by a finite system of ordinary differential equations. Given a local dynamical system, the ques- 
tion arises of finding all its gauge symmetries. If the equations of motion were variational, the 
Dirac-Bergmann algorithm [I], [2] would do the job. In fact, the Dirac-Bergmann algorithm does 
even more: it iteratively brings the variational equations of motion to a canonical form of the con- 
strained Hamiltonian dynamics with a complete set of constraints on the phase-space variables. In 
this form, the equations of motion are totally self-contained, having no other constraints that can 
be derived from any compatibility condition. It is the normal form of the variational dynamics, 
which is exploited in physics for many different purposes. Besides the other things, this canonical 
form allows one to classify the constraints by grouping them into the first and second classes, 
and it also results in identifying the gauge symmetry generators as the Hamiltonian vector fields 
associated to the first-class constraints. For the general (i.e., not necessarily variational) local 
dynamics, no canonical form has been yet identified which might be viewed as an equivalent of 
the Hamiltonian dynamics with a complete set of constraints. The correspondence between the 
complete set of first-class constraints, including the primary and secondary ones, and the gauge 
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symmetries is known as the Dirac conjecture. This conjecture is not always valid even for the 
regular variational equations, and the counterexamples are well known [2]. Beyond the scope of 
validity of the Dirac conjecture, no systematic procedure has been known for identifying the gauge 
symmetries even for the constrained Hamiltonian dynamics. 

In this paper, we work out an algorithm of bringing the general (not necessarily variational) 
dynamics to a certain normal form. In this form, the gauge symmetries and the complete set of 
constraints are explicitly identified. Since no natural pairing exists between the constraints and 
gauge symmetries unless the dynamics are Hamiltonian, our algorithm derives them independently: 
the constraints are found first, and then the symmetries are identified. Despite of the fact that 
Dirac's classification of constraints is conventionally defined in terms of the Poisson brackets, the 
notions of the first- and second-class constraints can be naturally extended to the non-Hamiltonian 
systems. The algorithm works equally well in the constrained Hamiltonian dynamics, finding out 
all the gauge symmetries, even though the Dirac conjecture does not hold for the system. So, 
the method may have some new impact on the well-studied area of the constrained Hamiltonian 
dynamics. 

Below we give some introductory comments on the background of the problem and outline the 
contents of the paper. 

Any system of ordinary differential equations can always be depressed to the first order by 
introducing new variables, so we start with a first-order autonomous system. Under certain 

n 

regularity conditional the first order ODE system can be always brought (by adding new auxiliary 
variables, if necessary) to the following form 

(1) x* = V*(x) + X a Z i a {x) , a = l,...,l, z = l,...,n, 

(2) T a (x) = 0, a = l,...,m, 

where x l , X a are the decision variables. Thus, in the regular local dynamics, it is always possible 
to have only two types of variables: (i) the phase space coordinates x l being subject to the normal 
differential equations (pQ) and the constraints fl2]), and (ii) the variables X a entering linearly and 

Prom the viewpoint of physics, the regularity means that the system has a definite number degrees of freedom, 
i.e., the equations admit the same number of independent Cauchy data in every domain of the configuration space. 
A more accurate formulation of the regularity conditions is given in Section 2 as well as the explanation how to 
bring the general regular equations to the normal form (Tlj), (|2"T). 
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without derivatives into the right hand sides of the differential equations for x's. We call the set 
of equations (CQ), (j2J) the primary normal form of local dynamics. Obviously, the compatibility 
conditions between the differential equations ([I]) and the algebraic constraints (j2J) can result in 
more algebraic equations, which are called the secondary constraints. The differential consequences 
of equations ([I]), ([2]) are considered in Section 3. 

Two particular classes of equations ([!]), ([2} were extensively studied by two different sciences: 
Dirac's analysis of the constrained Hamiltonian systems and optimal control theory. The first one 
deals with the case where the number of the constraints coincides with the number of A's, i.e., 
/ = m, the vector fields Z a and V are all Hamiltonian, and 



Here {■,■} is a non-degenerate Poisson bracket and T a are the algebraic constraints (T5]). Upon 
this identification, the equations of motion become variational, following from the least action 
principle for the functional 



where pidx 1 is a symplectic potential associated to the Poisson bracket. This can be understood 
as a conditional extremum for the standard Hamiltonian action subject to the constraints (J2]), 
with the variables A being the corresponding Lagrange multipliers. It is the variational dynamics 
that the Dirac-Bergmann algorithm is normally applied to. If, after applying the algorithm, the 
theory appears containing the first-class constraints, the solutions are not unique for the Cauchy 
problem and the action (HI) possesses gauge symmetries. 

Optimal control theory deals with another particular case of the above equations, which is 
opposite, in a sense, to the variational case: only the differential equations (0Q) are considered, 
with no constraints (T2]) imposed on x's. In optimal control [3], [4], the A's are called the control 
functions. If only equations ([1]) are considered, the Cauchy problem for x's will have a unique 
solution corresponding to every choice of the functions X a {t). The control functions remain un- 
restricted anyhow by equations ([1]), then the solutions remain ambiguous for x l unless all A's are 
determined by some extra requirement. The distinctions between Dirac's constrained dynamics 
and optimal control theory go far beyond imposing (or not imposing) constraints and supposing 
(or not supposing) the existence of a symplectic structure: they profess quite opposite concepts of 
the ambiguities in solutions of equations (CEJ). Optimal control theory considers the functions A as 



(3) 



K = {** , T a } 



V 1 = {x i , H } . 



(4) 
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describing a background for the dynamics, making it, in fact, non-autonomous. The objective is 
to minimize certain function(al)s of the solutions x % (t) (the cost functions) by varying the control 
functions A(t). Quite opposite, the Dirac constrained dynamics consider all the solutions for x(t) 
and X(t) equivalent to each other whenever they correspond to the same Cauchy data. In partic- 
ular, a function(al) is considered as physically observable, if it evolves in the same way on every 
trajectory from the equivalence class, i.e., it is out of control from the viewpoint of optimal con- 
trol theory. The gauge theory treats A as dynamical variables, not the functions controlled from 
outside. From this viewpoint, the controllable values are supposed to be unobservable because 
the autonomous physical models must be self-contained, leaving no room for intervention of any 
ultramundane force that can control/optimize what we observe. With this regard, the ambigu- 
ity brought to the solutions by the undetermined multipliers A(i) is to be factored out from the 
dynamics, first classically and then quantum-mechanically. A basic tool for such a factorization 
is a gauge symmetry relating equivalent trajectories. In particular, the physical observables are 
identified with the gauge invariant function(al)s. 

By a gauge symmetry of equations (JTJ, (j2J) we understand the infinitesimal transformations of 
the form 

v p+i 

(5) 5 £ x l = YR) y , 5 £ \ a = YK . 

n=0 n=0 

(n) 

The transformation parameters e are arbitrary functions of time, and e stands for the n-th 
order time derivative of e. The structure functions R\ n) and Ufa are allowed to depend on a finite 
number of variables x\ \ a , X a , X a , ... . The transformation flSJ) is supposed to leave equations (pQ), 
§2§ invariant in the sense that the variations 

(6) 5 £ (x l - V\x) - Z'Ux)\ a ) , S e T a {x) 

must be proportional to the equations of motion (CTJ, (j2J) for arbitrary e(t). In what follows we use 
the ordinary physical terminology and refer to the values that are proportional to the equations 
of motion as vanishing on-shell. 

Consider, for example, unconstrained equations (CQ) and suppose that the vector fields Z a span 
an integrable distribution of rank m, 

(7) [Z a ,Z p ] = UT J {x)ZT, 
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which is preserved by V, 

(8) [Z a ,V] = V£{x)Z p . 

Then there are exactly m linearly independent gauge transformations of the form ([5]). These read 

(9) 5 £ x l = Zle a , 5 £ X a = e a - (V£ + U^) . 

Notice that equations ([1]) are fully consistent, no mater what are the vector fields V and Z a 
involved in the right hand side. The involution relations ([7]), ([8]) are not the necessary conditions 
for consistency of the equations. If the distribution Z = span{Z a } is not integrable and/or not 
preserved by V, the explicit form of the gauge transformation ([5]) has remained unknown yet, and 
we are going to present it in Section 4. 

Given the distinctions in the main concepts/objectives of the gauge theory and optimal control, 
it is no surprise that the gauge symmetry remained unstudied for the general dynamical equations 
(CQ), (T5]) in the context of optimal control. Even though a lot of interesting non- variational gauge 
models arise in the context of modern higher energy physics, e.g. self-dual Yang-Mills fields, 
higher-spin gauge theories, Sieberg-Witten and Donaldson-Uhlenbeck-Yau equations, etc, the issue 
of extending the Dirac-Bergmann method to the general local dynamics has been shelved yet, 
because no ways were seen to quantize the non- variational dynamics. Recently, we have proposed 
the BRST quantization methods applicable to arbitrary non- variational dynamics [5] (see also 
[6] for another version of this method). Our method implies the characteristic distribution to 
be on-shell involutive and tangent to the constraint surface. As it has been already explained, 
this is not the most general case. In Section 5, we explain that any local gauge dynamics can 
be eventually reformulated in the involutive form, which is physically equivalent to the original 
one. The involutive normal form of dynamics, resulting from the extension of the Dirac-Bergman 
algorithm to non- variational dynamics, thus becomes a starting point for applying the deformation 
quantization techniques for the general dynamical systems. 

2. Regularity conditions and the primary normal form of local dynamics 

Let us start with explanation how to bring a local dynamical system to the normal form (TjQ), 
(T5]) . Any finite system of ordinary differential equations can always be depressed to the first order 
by introducing new variables. By further adding new variables, if necessary, the first-order ODEs 
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can be equivalently rewritten in the form of inhomogeneous Pfafnan equations 

(10) Oj^x 1 = Vj(x), J=1,...,N, z = l,...,n, 
subject to constraints 

(11) T a (x) = 0, a = l,...,m. 
For example, given second-order equations 

(12) f{y,y,y) = o, 

by adding new variables v and w, these can be equivalently rewritten in the form ( jTUl) . ( lTTj) as 

(13) v = w, y = v, f(y,v,w) = 0. 

Sometimes, to take into account the geometric properties of the original equations (TT21) . it might 
be more appropriate to introduce the new variables v , absorbing the derivatives of x, as linear 
inhomogeneous functions of x, v — e(x)x + T(x). This will result, however, in the equations of the 
form (fTU| [TTI) . anyway. We call the primary constraints or just constraints both equations (TTTT) 
and the functions T a (x). 

To further proceed with the general equations of motion in the constrained Pfaffian form (fTUl) . 
(TTTj) . we need to impose appropriate regularity conditions. The phase space of the dynamical 
system is supposed to be an open domain U in R n with linear coordinates {x 1 }. All the functions 
of x % entering equations (ITUl) . ffTTj) are supposed to be analytic. E 

Let us introduce some convenient notions for describing the regularity of the equations. Let M 
be a matrix whose elements are analytic functions on R™. Given a subset [/ci", one can define 
a sequence of embedded subspaces 

(14) U= U° D U l M d U 2 M D ■■■ D E/S = 
where 

(15) U r M = {xeU\ rankM(x) > r} 



2 The conditions we adopt here are not the weakest possible ones. For example, one can extend the consideration 
from W 1 to an analytic or smooth manifold M, with constraints {T a (x)} being a section of a vector bundle over 
M. In this paper, however, we are going to avoid the technicalities related to less restrictive regularity conditions, 
for the sake of a more clear presentation of the algorithm as such. 
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is the set of all points in which the rank of the matrix M is not less than r, and m is the maximal 
value of rankM on U. Then V r M is an open everywhere dense subset of U. A minimal subset 
U' M C U will be called a a regular part of U regarding M. It will be denoted just by U' whenever 
the corresponding matrix- valued function M is obvious from the context. The complementary 
subset to the regular part, U\U', will be called an abnormal part of U. 
Denote by £ the common zero locus of the constraints (iTTj) . i.e., 

(16) £ = {x E U\ T a (x) = 0,a = l,...,m}. 

The primary constraints are assumed to be consistent as algebraic equations per se. If the con- 
straints contradict each other, the entire system of the dynamical equations flTJ), (J2J) does not have 
any solution, and this can hardly be considered as regular dynamics. The same consistency of 
the constraints by themselves will always be assumed for the secondary constraints derived in the 
next section. 

Let U' be the regular part of U regarding the Jacobi matrix of the constraints J = (diT a ). The 
constraints {T a } are called regular if 

(17) E' = [/'n£/0. 

In this case, £' is a smooth submanifold in U' with dimS' = n — rankj|[/'- 

Denote by E" the regular part of £' regarding the matrix = (9ji) of the Pfaffian forms in 
( flOl) . In accordance with the definitions above, E" C U = (U')q, where r = rankB(p) for p e 
Equations (flOl) are linear in velocities x l , with the coefficient matrix (6ji) having rank greater 
than or equal to r on U. Restricting the dynamics to U, we can solve (FlOl) with respect to the 
velocities as 

(18) x i = V\x) + \ a Z l a (x) , a = l,...,Z, l = n -r. 

Here X a (t) are arbitrary functions of time, the vector field V is any solution to the inhomogeneous 
linear equations 

(19) (djiV'-v^y^o, 

and the vector fields {Z a } span the space of solutions to the corresponding linear system 

(20) 0j;ZV = O. 
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Upon restriction to the domain U, the system of equations ([TBI) and (fill) is completely equivalent 
to the original equations (flOl) and (lllj) . We refer to equations (I18j) . ( TTTj) as a primary normal 
form of local dynamics. It is the form ([1]), (T5]) announced in Introduction. The vector field V 
will be called a primary drift, and the vector distribution Z = span{Z a } will be called a primary 
characteristic distribution. Notice that the distribution Z is not required to have a constant rank 
on U nor is it supposed to be involutive. 

The primary normal form of equations (00), (J2J) implies certain equivalence relations among their 
ingredients. Besides the non-degenerate changes of the decision variables 

(21) x l ^ x H = f{x) , X a ^\ ,a = F°(x)\> 3 + G a (x), 
the system is invariant with respect to the following redefinitions: 

(22) T a ^T' a = S b a T bl V^V' = V + T a W a , Z a ^ Z' a = Z a + T a X a a , 

where {S h a (x)) is an arbitrary non-degenerate matrix on U, and W a and 1° are arbitrary vector 
fields. So, neither the primary drift V, nor the characteristic distribution Z is defined uniquely by 
the inhomogeneous Pfaffian system (fT0|) subject to the constraints (111!) . We say that two dynamical 
systems on U are equivalent if their primary normal forms (pQ), ([2]) are related by transformations 
( 12T1) . ( |22i) . In what follows, we will deal with those dynamical systems which admit at least one 
analytic representative in each equivalence class, and only the analytic representations will be 
considered. 

A function (or vector field) F is called trivial if it vanishes on With account of regularity 
we have 

(23) F\y,' = <£> F — T a F a 

for some smooth functions (or vector fields) F a . From this viewpoint, the automorphisms ( J2T1) . 
(]22|) imply that two characteristic distributions are equivalent whenever their difference is trivial, 
and two drifts are equivalent whenever they coincide modulo the characteristic distribution and a 
trivial vector field. 

As it has been already mentioned in Introduction, the general equations of motion, being brought 
to the primary normal form (CQ), (E]), can have nontrivial compatibility conditions resulting in 
additional constraints on x's and A's. The secondary constraints on x's can be derived directly in 
the inhomogeneous Pfaffian form ( FlOl) . (TTTT) . as it was shown by X. Gracia and J.M. Pons [9]. In 
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the next section, we consider the procedure of deriving the secondary constraints and fixing the 
undetermined multipliers A by making use of the primary normal form of the equations of motion. 
This procedure allows us to extend the Dirac classification of constraints to the general equations 
(CQ), (J2J), without appealing to any Poisson structure. 

3. Consistency conditions and secondary constraints 

The consistency conditions of the primary normal equations originate from the requirement 
that the integral trajectories of the flow (DO), with fixed functions A(i), should be confined at the 
regular part of the constraint surface (TTTF) . For the regular constraints this means that the time 
derivatives of all the constraints must vanish on the integral trajectories whenever they intersect 
E': 

(24) T a | on _ shell = (VT a + \ a Z a T a ) \v = . 

The expression in brackets being a trivial function in the sense of ( 123]) . we arrive at the following 
linear inhomogeneous equations for A: 

(25) VT a + \ a Z a T a = F b a T b . 

Let E" be the regular part of E' regarding the matrix M = (Z a Tb), and rankM(p) = s for p e E". 
Then, according to the definition (fT5|) . E" C (U') S M . 

To further proceed with solving equations (T251) . we decompose the constraints, characteristic 
distribution and undetermined multipliers in the following way: 

(26) T a = (T A ,T Sl ), Z a = (Z A ,Z ai ), X a = (X A , A" 1 ) , 

(i) 

where A = l,...,s and D = [Z aTb) is a maximal non-degenerate minor of M. Denote U the 
regular part of (U') S M regarding D. Using equations (125]) with a = A, we can express \ A as linear 
inhomogeneous functions of A Ql : 

(27) X A = -D AB (VT B + X ai Z ai T B ) , 

with the matrix (D AB ) being inverse to D. After substitution of the determined multipliers ( 1271) 
into the remaining equations (125|) . we get 

(28) T}+\^B ai - ai =0, 
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where 

(29) 2* = VT- ai - D AB (VT B )Z A T- ai , B ai&1 = Z ai T- ai - D AB (Z ai T B )Z A T- ai . 

Obviously, the matrix B = (B ai5il ) must vanish on since otherwise one could determine more 

multipliers A, that would contradict maximality of the non-degenerate minor D. Being a trivial 

(i) 

analytic function in JJ, the matrix B has the form B aiai = W£ a T c . Equations (|28|) should be 
satisfied together with the equations of primary constraints, hence the trivial terms containing A's 
can be omitted and we obtain the new constraints 

(30) Tl(x)=0. 

So, we see that the conservation laws (124p for the primary constraints are equivalent to equations 
(12"?]) determining a part of A's as specific functions of x's, and the A-independent relations (!30|) . 
Equations (1301) as well as the functions T£ themselves are called the secondary constraints of 
the first stage. Notice that the secondary constraints are defined modulo trivial contributions 
proportional to the primary constraints. Also notice that the functions T x (x) are not necessarily 
independent. 

Substituting the fixed multipliers fTSTj) into the original equations of motion (pQ), (T5]) and adding 
the secondary constraints, we arrive at the following set of equations: 

(31) x l =V l (x)+X^ Z^ix), T ai (x) = 0, 
where the values 

(32) V= V - D AB (VT B )Z A , Z ai = Z ai - D AB (Z ai T B )Z A , T ai = (T a , T^) . 

are called, correspondingly, the drift, characteristic distribution, and constraints of the first stage. 

(i) 

By construction, all these objects are well defined on the open everywhere dense domain JJC U . 
Let us write down the following obvious relations: 

(33) V T A = , V T- ai = 2* , Z ai T A = , Z ai T ai = W c aiai T c . 

Along with the condition of invertibility of the matrix D = (Z A T B ) these relations give a tip 
for a simple geometric interpretation of the transition from the original equations ([IT), (l2|) to 
the equivalent first-stage equations (13TT) . The original surface of primary constraints ( TT61) can be 
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represented (at least locally) as a transverse intersection E = E 11 n S ± of two surfaces E and E x . 
These are given by 

(34) E 11 = {xeU | T ai {x) = 0}, S X = {x€c7 | T A (x) = 0}. 

In its turn, the original characteristic distribution Z = span{Z a } is decomposed into the sum 
Z = Z n © Z ± of two sub-distributions 

II (1) _L 

(35) Z =sp&n{Z ai }, Z = sp&n{Z A } , 

which are called, correspondingly, the tangential and transverse. The multipliers X A , being related 
to the transverse sub-distribution, are fixed by the conservation law of the primary constraints, 
whereas the multipliers A" 1 corresponding to Z^ still remain undetermined. 

Geometrically, relations (|33|) mean that the distribution is tangent to the regular part of 
the primary constraint surface E, whereas the complementary distribution Z^ is transverse. The 
tangential distribution Z^ C Z is invariantly defined by the property Z^ |s C TE. The dimension 
of the transverse distribution, being complementary to the dimension of z\ coincides with the 
number of the multipliers (1271) determined from the conservation requirement for the primary 
constraints. The primary characteristic distribution Z is tangent to the surface E 11 , being zero 

locus of the primary constraints not involved in determining of the multipliers. Also notice that 

(i) 

the first-stage drift V is tangent to the primary constraint surface E on the zero locus of the 
secondary constraints. 

Not only do the first-stage equations (1311 describe the same dynamics as ([T]), (j2J) (the zero-stage 
equations), but they appear identical to equations ([ID, (j2J) in form. Therefore, we can apply to 
( 13~TT) all the above reasonings concerning the preservation of the constraints in time. Supposing 
the first-stage constraints to be regular, we check whether the conservation condition 

(36) A. $ I = o 

v ' dt 1 a i Ion-shell 

implies new secondary constraints, or it further restricts the undetermined multipliers A. The new 
constraints, if any, are called the secondary constraints of the second stage and denoted by T| . 

(2) 

We suppose that the complete set of the second-stage constraints T a2 = (T a ,T± i: T? 2 ) is regular. 
After exclusion of the determined multipliers from the differential equations, we get the drift and 
the characteristic distribution of the second stage. If there are nontrivial secondary constraints 
among the constraints of the second stage they must conserve. So, we further proceed with 
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deriving consequences from the conservation of the second-stage constraints. The algorithm will 
continue working, defining some multipliers and redistributing the corresponding vector fields from 
the tangential distribution of the previous stage to the transverse distribution of the next stage, 
redefining the drift, and bringing the constraints of the next stage. Since the overall number of 
the independent regular constraints at any stage cannot exceed the dimension of the phase space, 
the algorithm has to stabilize after a finite number of steps. The stabilization is achieved as soon 
as new secondary constraints stop to appear from the consistency conditions. If the iterative 
procedure terminates at the k-th stage, the equations of motion take the following form: 

(37) x i =V i (x)+X^Z i ak (x), T ak =0. 

These equations are defined in an open everywhere dense domain JJC U, and the complete set of 
constraints 

(38) T ak = (T a ,Tl,---,Tl) 

includes the primary and secondary constraints of all stages. 
By construction, the following relations take place: 

(fcj(fc) b (k) (k) (k) b w 

(39) VTa k — F ak Tb k , Z a k T a k = F a h k a k T b k 

for some functions and F^ a . We call (EJT1) a complete normal form of a local dynamical 
system. It is the form that ensures full consistency of the dynamics as it does not have any further 
compatibility conditions and consequences. 

Several remarks can be made about the algorithm above: 

• The (n+l)-st step of the algorithm becomes possible whenever the constraints of n-th stage 
are regular. Actually, it is the constraint surface that has an invariant geometric meaning, 
not the constraints as functions. The same constraint surface can be defined by different 
T's, and all such constraints are considered to be equivalent. If the equivalence class 
includes a regular representative, this representative is picked up as the set of constraints 
of the n-th stage. As is seen from the iterative procedure, the algorithm is not sensitive to 
the specific choice of regular representative. 

• The (n + l)-st iteration involves the ambiguity concerning the choice of a maximal non- 
degenerate minor D n of the matrix M n = [Z an T a „)- The different choices can result 

(n+l) 

in different domains of definition JJ for the system of the (n + l)-st stage. Notice that 
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(i+i) 

the different domains U s, being open everywhere dense subsets in U, should coincide 
modulo a set of measure zero. In the intersection of these almost coinciding domains of 
definition, the equations of motion are not sensitive to a specific choice of a maximal minor. 
• When the compatibility conditions (125]) are considered as a system of linear equations for 
A's, the surface S" was defined as a subset of £', where rankM = s is maximal. This 
requirement of maximality can be relaxed, in principle. One can consider a subset ScE', 
where rankM = s' is less than s. According to the terminology of Section 2, S belongs to 
the abnormal part of the constraint surface £' regarding M. In some cases, the abnormal 
set S C U' can be a smooth submanifold defined by a set of regular constraints T s = 0. 
At every given stage, the dimension of the abnormal part of the constraint surface is less 
than the dimension of the regular part, so that dimS 1 < dimS. Being restricted to S, the 
dynamics would be regular in the same sense as with the restriction to the regular part 
of the constraint surface. With the restriction of the dynamics to the abnormal set S, a 
lesser number s' of the multipliers A are determined, but the greater number of the con- 
straints appear at this stage. Thus, at any stage of the iterative procedure, there may be 
an ambiguity in deciding between the regular or abnormal parts of the constraint surface. 
The dynamics can be regular and consistent for either of these options. It is amply clear 
that the dimension of the complete constraint surface (J38l) and the number of undeter- 
mined multipliers in fl37j) can depend on the choice made at certain stage (the algorithm 
"bifurcates"). So, after applying the algorithm, which works equally well with regular and 
abnormal parts of the constraint surface, one can learn that the original dynamical system 
(CQ), (J2J) has contained in itself several different dynamics (1371) . with different constraints, 
drifts and numbers of undetermined multipliers. These dynamics never communicate and 
should be interpreted as different physical systems, even though the difference reveals itself 
only after applying the algorithm, being not explicitly visible in the primary normal form 
of the dynamics. 

Let us comment on the geometry of the dynamical equations in the complete normal form (|37|) . 
The primary characteristic distribution Z = span{Z a } is eventually decomposed into the direct 
sum 



(40) 



Z = Z± © Z\ 
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of the transverse and tangential distributions with respect to the complete constraint surface 

(41) Z = {xEU | T ak =0}. 

(fe) 

As is seen from the tangential distribution = span{Z Qfc } preserves the complete con- 
straints surface ( I4TI) . and the multipliers A Qfc related to the basis vector fields of Z\\ remain arbi- 
trary functions of time in equations (137|) . The distribution Zx ~ Z/Z\\, being a complement to 
J?y, is transverse to the complete constraint surface E: If {Z^} is a basis in Zx, then 

(42) rank(Z^ T J = dim Zx. 

All the multipliers X A corresponding to Z± are defined by the conservation conditions of the 
primary or secondary constraints, at one or another stage of the algorithm. 

By making use of the equivalence transformations ( 1221) . the complete set of constraints can be 
rearranged into the union of two subsets 

(43) T ak = in , r£) 
such that 

(44) Zfl\ t = MZeZ and detZ>| s ^0, 

where D = (Z^T^). Denoting by Sy and Ex the zero loci of the constraints T 11 and T x , respec- 
tively, we see that the complete constraint surface can be represented as the intersection of two 
surfaces: 

(45) £ = En n S± . 

The primary characteristic distribution Z is tangent to Ey on E and transverse to Ex. Notice 
that codimEx = dimZ±, while the dimension of the tangential distribution Z\\ does not correlate 

with the dimension of Ey in general. 

(fe) 

The complete drift V = V defined by f l37|) is constructed step by step from the primary drift V 
by adding terms proportional to the elements of the transverse distribution. When the algorithm 
has terminated at the k-th stage, the drift reads 



(46) V = V - {VT^)(D' l ) AB Z 



B 



Notice that V belongs to the equivalence class of the primary drift V in the sense of the equivalence 
relations (1221). 
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Let us elaborate on the correspondence between the basic notions of the constrained Hamiltonian 
mechanics and their counterparts in the general local dynamics. The equations ( 1371) can be 
viewed as an equivalent of the Dirac constrained dynamics with the complete set of primary and 
secondary constraints ( 1381) . The "tangential" and "transverse" constraints T 1 ' and correspond, 
respectively, to the complete sets of the first- and second-class constraints in Dirac's classification of 
the Hamiltonian constraints. The tangential and transverse primary characteristic distributions 
Z\\ and Z± correspond to the distributions generated by the Hamiltonian vector fields of the 
primary first- and second-class constraints, respectively. One may also wonder about an analog 
for the Hamiltonian vector fields associated with the first-class constraints. These are known 
to generate the whole set of gauge symmetries provided that the Hamiltonian system obeys the 
Dirac conjecture. From the next section we will learn that the vector fields in question form the 
distribution 

(47) Z v = Z l{ U[Z h Z l{ ]U[Z h V]U--- , 

where the dots stand for the higher iterated commutators of Z\\ and V. Whether the Dirac con- 
jecture is true or not, it is the distribution Zy that generates the gauge symmetry transformations 
(JSJ) of the dynamics. 

4. Gauge symmetries 

In the previous section, we have seen that the algorithm of stabilization of the primary con- 
straints brings the original equations of motion to the complete normal form fl37j) . (1391) . In this 
section, we find all the gauge symmetry transformations for these equations. 

To ease the notation, we omit all the sub- and superscripts referring to the final stage of the 
iterative procedure from the previous section. After this omission, equations ( 1371) read 

(48) x l = V i (x)+\ a Z t a (x), T a (x) = 0, 

where the distribution £ = span{ V, Z a } is assumed to be tangent to the complete constraint 
surface E = {x 6 U\ T a (x) = 0}. 

The infinitesimal gauge symmetry transformations (jHJ) are sought for equations (l4"8l) in the form: 

p p+i 

(49) 8j = £ 4-n * > ^ = E U P^-n ( Z > 

n=0 n=0 
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. .. (A;) 

where i?'s and £/'s are some functions of x, A, A, ... up to some finite order A, and e(t) are the 
transformation parameters, being arbitrary functions of t. Notice that the integer p, called the 
order of gauge transformation, is not fixed a priori. Moreover, it can always be risen by the trivial 
replacement e i— > e of the gauge transformation parameter. Conversely, if R p = and U p+ \ = 0, 
then the reverse change of variables e i— > e depresses the highest order of the derivative of e in 
( 1491) . When the functions R p and U p+ i do not vanish simultaneously^, the highest order of the 
derivative cannot be depressed and we refer to corresponding transformation as undepressible. 

Let us first consider the issue of gauge symmetry for the dynamical system (148j) without con- 
straints T a . The gauge invariance of the dynamics with respect to the infinitesimal transformations 
l9|) means that the gauge variation of (HHl) should vanish on their solutions, i.e., 



(50) 6 £ (x l -V*-\ a Z l a )\ on _ shcU = 0. 

Excluding the velocity x % in (1501) with the help of ( HHj) . we arrive at the following recurrent relations 
for the structure functions R and U of the gauge transformation (149]) : 

(51) Rl = W* , 

(52) i^+i = DRi n + w Ui j n = 0,...,p-l, 

(53) r; +1 = dr; + w; +1 = o . 

Here we have introduced the following abbreviations: 

(54) Wi = UZZ* a , D = -d-[V + X a Z a ,-], 

and d is understood as the time derivative acting only on A's, i.e., 

• <9 •■ d 

(55) d = \ a — + \ a ^ + ... . 

Equation (|5~1|) implies that the coefficient at the highest derivative of e in 5 £ x l is given by a linear 
combination of vector fields from the characteristic distribution Z. Relations (1521) enable one to 



3 Below, we will prove the following fact: For an undepressible gauge transformation all the vector fields Rk, 
k = 0, 1, are linear independent. This means that the gauge variation of x's contains actually all the derivatives 
£,i,..., e , without gaps. 
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express all the structure functions R in terms of the functions U defining the gauge transformations 
of the undetermined multipliers A. Namely, by solving relations ( |52i) . we find 



In particular, equation (1531) takes the form 

p+i 

(57) Rl +1 = J2D m W; +1 _ m = 0. 

m=0 

The last equation on U's is the only nontrivial condition to satisfy. Its solutions form a linear space, 
which dimension is to be computed, depending on p. To this end, certain regularity conditions 
must be imposed on the equations of motion and their domain of definition. Also, we will use 
some properties of and notions about the distributions, which are briefly listed below. 

So far, the components of the vector fields V, Z a were supposed to be analytic functions defined 
in some domain of linear space. In the following, to avoid specifying each time the definitional 
domain for the gauge transformations, we relax the condition of analiticity and utilize the algebraic 
setting usually adopted in control theory [3] for the dynamical equations jT]). Namely, V and Z a 
are allowed to be meromorphic vector fields. As the scar set of a meromorphic function is of 
measure zero, our gauge transformations will be well defined in an open everywhere dense domain 
U of the original phase space. It is U that is to be considered as the definitional domain of the 
gauge dynamics. 

Let F be the field of meromorphic functions and W be the space of meromorphic vector fields 
on U C M n . The space W has both the structure of a real, infinite-dimensional Lie algebra with 
respect to the commutator of vector fields and the structure of n-dimensional vector space over 
F. We will refer to F-linear subspaces of W as distributions. Generally, an arbitrary distribution 
P C W is not closed with respect to the commutator, i.e., the commutator of two vector fields 
from P may not belong to P. If [P, P] C P, then the distribution P is said to be involutive. 
The Lie closure of P is defined as the minimal involutive distribution P C W that contains P. 
Uniqueness of P follows immediately from the fact that the set of the involutive distributions is 
closed with respect to intersecting of linear spaces. Furthermore, one can inductively see that the 
distribution P is generated by the vector fields 



■n 



(56) 




(58) 



[• ■ ■ [ui,u 2 ],w 3 ], • ■ -],Uk] \/u k e P , VfceN. 
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For an involutive distribution P, P — P. 

By taking the multiple commutators ( l58i) . one can filter the Lie closure P by the sequence of 
sub-distributions 

(59) Pi = span^-j [• • • [m, w 2 ], u 3 ], ■ ■ ■ ], u k ] \ Wu { e P , k = l,...,l } 
such that 

(60) P = P m ^ Pm-l 3 P m - 2 D ■•OP 1 = P. 

This filtration is known as the Lie flag of the distribution P. Clearly, the minimal integer m < n — 
dim P involved in (160]) is an invariant of the distribution P together with the numbers df, = dim P^. 
The integer m is usually referred to as the depth of the distribution P and the sequence of integers 
(di, d 2 , d m ) is called the growth vector of P. 

To check the solvability of equations flSTj) we need some extension of the field F. Along with the 
coordinates x % in W 1 we introduce the infinite set of variables {A^}, a = 1, ...,m, fceN. Denote 
by T the field of meromorphic functions of x % and a finite number of the variables A£. Replacing 
in the definition of W the field F by its extension J 7 , we get the n-dimensional vector space 

The JF-Iinear subspaces of the vector space W will be called A-distributions. 
Let us turn JF into a differential field by setting 

dz* = , d\% = X a k+1 . 

Clearly, this definition just mimics the definition of the time derivative (1531) if one set Aq = X a . 
The action of the operator d can be further extended to the A-distributions by the rule 

d(aw) = (da)w + adw Va G T , Vw G W , 

With all the definitions above, we can start studying the gauge symmetries (j49p of the system 
( 1481) . Every dynamical system (HHl) defines and is defined by the distribution £ generated by the 
vector fields V and Z a . Associated to £ is a gauge distribution Z v . The latter is defined as the 
limit of a filtration 

(62) C C ■ ■ ■ C = Z v CW 
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of involutive A-distributions Z v given by 



(63) 2*=\J [V,...,[V,Z]. 

m=0 V v ' 

m 

It is clear that Z v = Z and B = Zy U V. 
Define an R-linear operator D : W — > W by 

(64) Dw = -dw - [V, w] - \ a [Z a , w] VwGW. 

Although the operator D is not ^-"-linear, it satisfies the following analog of the Leibnitz rule: 

(65) D(aw) = (Da)w + aDw Va G T , VweW, 
where 

(66) Da = -da -Va- \ a Z a a . 

Since D is a differentiation of the field J 7 , one can thought of W as a differential ^-"-module. 

Given an element w G Z, consider the sequence of elements D k w G Zy, k G N. As Zy is of 
finite dimension, there is p G N and a set of functions ai, - • • , a p G such that 

(67) £> p u> = ai/^w; + a 2 D p - 2 w + ■ ■ ■ + a p w . 



Proposition 4.1. Lei w £ Z satisfy equation ft67\ ), then there exist a sequence of elements 
wi, ...,w p G Z such that 

(68) D p w + D^wx + D p ~ 2 w 2 + --- + w p = 

and Wk = bkW for some bk G T . 

The proposition is proved by induction, making use of the Leibnitz rule ( 1651) . 
Corollary 1. For any vector field Wq from the characteristic distribution Z = span{Z Q } 7 equation 



(51) has a solution for some p. In other words, every basis vector field Z a generates a gauge 



transformation, so that the total number of independent gauge parameters e coincides with dim Z. 



An element w G Z is said to have degree not higher than p, if it satisfies equation ( 1681) with 
some (not necessarily linear independent) Wk G Z. It is clear that the elements of degree not 
higher than p form a A-distribution Z p C W. We have the finite filtration 

Z = Z N D Z N ^ D • - • D Zx D . 
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The numbers 

5 P = dim Z p — dim Z v _\ 

are called the indices of the characteristic distribution Z. 

By making use of the Euclidean metric in MJ 1 , we can split the imbedding Z v _\ C Z p as 
Z p = Z p ^i © Z p _ v As a result the characteristic distribution is decomposed into the direct sum 

(69) Z = Z 1 @Z 2 @---@Z N , Z P ~Z P /Z P _ U dimZ p = 5 p . 
In a driftless theory V = and 5± > 1 as 

-D(AqZ q ) = A^-Zq, 

and equation (1681) is satisfied with p = 1. The corresponding gauge transformation (149j) is the 
time reparametrization. 

We use the decomposition ( 1691) to construct a basis of undepressible gauge transformations 
for the equations of motion (HHj) . If {Z ap } 5 ^ p=l is a basis in the A-distribution Z p , then the 
undepressible gauge transformations read 

(70) <w = £ ^;_„_ 1)Qp ( ? ap , u« = e ui_ n)ap f <* , 

n=0 n=0 

where 

n 

f 71 ") /?* = \ D m W i W* = JJ a 7 l JJ a = ft a 

^ I (»)% (n-m)<*P ' {n) a P (n) «P a ' (0)"p «P ' 

m=0 

Proposition 4.2. If w & Z p , then the vector fields w, Dw, D 2 w, D p ~ 1 w are linearly inde- 
pendent. 

This is true, otherwise it would be to e Z ? with q <p, because of Proposition I4.ll 



Corollary 2. TTie vector fields R\ n)ap entering $7u\) are linear independent, so that 5 £ x % involves 
all the successive derivatives of the gauge parameter up to the order k. 

Let us denote 

Z D = span^{D m n | Vu G Z, m E N} . 

By construction, the A-distribution Zp is invariant under the action of D, i.e., DZp C i^o- It 
turns out that the distribution Zp is actually involutive and coincides with Zy. 
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Proposition 4.3. Zb = Zy . 
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Corollary 3. The gauge variations (TUt) of x's are spanned by Z 



>v- 



Proof. Since Zd C Zy, it remains to prove the converse inclusion. Let us first consider the special 
case where V = 0. For a driftless system Zy = Z and we must show that Zjj contains the Lie 
closure of Z. 

(k) 

Increasing, if necessary, the order of the derivatives of the gauge parameter by changing e ^ e , 



the gauge transformations (170]) can be brought to the form 

(72) Sex* =e "tfjx) + ■ ■ • , 5 £ X a =T <* + ••-, 

where the dots stand for the terms with the lower order derivatives of the gauge parameter. As 
the transformations (1721) exhaust all the gauge symmetries of the equations, they have to form an 
on-shell closed gauge algebra with respect to the commutator of infinitesimal transformations: 

(73) [6 S1 , 8 £2 ] | on-shell = 4 3 ) 

where 

(74) 4 = £^Wl, m (*,A fc ) ( ?V, 

n,m 

with /'s being the structure functions of the gauge algebra. On the other hand, 

(75) [6 £l ,6 £2 }x'= { e )a 1 { e ) ^Z a ,Z^ + -- - , 

(fc) 

where the dots stand for the other bilinear combinations of the derivatives e " 2 . At every fixed 

(*0 

instant of time, the derivatives of the parameters e f 2 can take on arbitrary predetermined 
values. Then, comparing the right and left hand sides of (173|) . we conclude that all commutators 
[Z a , Zp] are given by linear combinations of elements from Z^. A similar analysis for the successive 
commutators of the gauge transformations 

shows that all multiple commutators 

are also included into Zp. In other words, we see that Z C Zp. 
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The general case of a non-vanishing drift, V 7^ 0, can be formally reduced to the previous one 
by the following trick. Let us associate to equations (148p another dynamical system 

(76) x l = eV\x) + X a Z a (x) , 



where e is a new undetermined multiplier. The system ( |76|) is driftless, and hence reparametriza- 
tion invariant. The characteristic distribution of (1761) is generated by the vector fields V and 
Z a . The basis of infinitesimal gauge transformations for this new system can be chosen in the 
following way. First of all, the system is invariant under reparametrizations. The corresponding 
gauge transformation reads 

(77) 8 e x* = ex 1 = e(eV i + \ a Z a ) , 5 t X a = eX a + eX a , 5 e e = ee + ee . 

As is seen the generator of this transformation involves the vector field V and acts nontrivially 
on the new multiplier e. Considering e as a fixed function of time, we then define the gauge 
transformations of the form 



(78) 5 £ x l =e a Z i a + ■ ■ ■ , 5 £ X a = T ] " + •••, 5 £ e = 



The existence of these transformations easily follows from Proposition 14.11 According to the 
general formulas flTOj) . fTTTl) the expansion coefficients of 5 e x l in the time derivatives of e's are 
given by linear combinations of the vector fields D™Z a , where 

D e = -d-e[V, - }-X a [Z a , •]. 

Since the leading terms in the variations 5 e x l and 8 £ x l span the entire characteristic distribution 
S = span-jV, Z a }, the transformations (1771) . (ITS]) exhaust all the gauge symmetries of (1761) . 

We claim that the ^-transformations (J7H1) constitute an ideal of the algebra of all gauge trans- 
formation (EZD, (J7HD, i.e., 

(79) [6 £1 ,6 £2 }\ 

on— shell ^£3 > on— shell 

To prove this statement it is enough to apply the commutators in the l.h.s. of the relations above 
to e. We have 

(80) [6 ei ,6 ea ]e = 0, [<f e ,<y ei ]e = 0. 

Thus the reparametrization transform (177jl does not contribute to the r.h.s. of ( 1791) . 



As before, we have risen here the order of the gauge transformations to a certain uniform value 
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Arguments similar to those we have used in the driftless case allow one to prove that all the 
multiple commutators of the vector fields V and Z a are contained in Zp. Consider, for example, 
the following commutator of the gauge transformations: 

(81) [5 e , 6 £ }x l = e e a e[V, Z a ]x l + ■■■ . 

Because of relations (1791) the right hand side of (IHTj) must be a gauge transformation of the form 
flTJJ. We have 

(82) [S e , 5 £ }x l = 5 £ ,x l , s"i = e £ a f2 n {x, e s , \ k ) + ■ ■ • , 

where the dots stand for the other combinations of the time derivatives of e and e a . Comparing 
the coefficients at e e in ( 1HT1) and ( 1821) . we see that every commutator e[V, Z a ] is given by a linear 
combination of D™Z a . Setting e = 1, we conclude that [V, Z a ] e Z D . 

In a similar manner, one can see that all the higher iterated commutators of Z a and V belong 
to Z D . 

□ 

By construction, the structure functions (|7I]1 are meromorphic functions of x % and A a ■ So, the 
trajectories (x l (t), \ a (t)) can exist such that the transformations ( 1701) are ill defined. It is easy to 
see, however, that the functions (I7TT) can be made real analytic by an appropriate change of the 
gauge parameters: 

(83) e a -»■ e a = U%e p , 

with det(C/g) being a nonzero element of T . The real analytic transformations are well defined 
on all the trajectories. 

Example. Consider the driftless system 



(84) 



x* = rzi 
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associated to the following characteristic distribution in M 10 : 



1 ® 4- x 2 ® + x 3 ® + x 4 ® 
1 dx 1 dx 8 dx 9 dx w 

d o d 4 d 
Z 2 = ^+x 3 — + x 4 



(85) 

^3 



dx 2 dx 6 dx 7 
Z, = ^ + x 4d 



Z 



dx 3 dx 5 ' 
d 



i 



dx A 



A straightforward computation shows that the vector fields {Z a , [Zp, Z 7 ]} span the entire tangent 
space of M 10 and 

[[Z a ,Z p ],Z 7 ] = 0. 

Introduce the vector field X n = X^Z a and the differential D = —d — [X , ■ ]. One can readily 
check that 

DX + X 1 = , 

D 2 X 1 + 2DX 2 + X 3 = , 

(86) 

D 3 X 2 + 3D 2 X 3 + 3DX 4 + X 5 = , 

L> 4 X 3 + 4L> 3 X 4 + QD 2 X 5 + ADX 6 + X 7 = . 



These equalities are the particular cases of the general identity 



(87) ^CtD m - fc X m+ ^ 1 = 0. 

k=0 



The vector fields {X , X l5 X 2 , X 3 } form another basis in the A-distribution span^ZoJ. Com- 
paring relations (1861) with (1701) . (17B and (1571) . we get the following gauge transformations for the 
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differential equations (1841 : 

5 £ x = ExXq 

+ i 2 X 1 + e 2 (DX 1 + 2X 2 ) 

+ e 3 X 2 + e 3 (DX 2 + 3X 3 ) + e 3 {D 2 X 2 + 3DX 3 + 3X 4 ) 

+ e 4 X 3 + e 4 {DX 3 + 4X 4 ) + e 4 {D 2 X 3 + ADX A + 6X5) + e 4 (D 3 X 3 + AD 2 X 4 + 6DX 6 + 4X 6 ) , 

S £ \ = e 1 \ + £ 1 \ 1 

+ e 2 Ai + 2e 2 \ 2 + e 2 X 3 

+ e 3-^2 + 3£3A3 + 3e 3 X 4 + E 3 \§ 

+ "e\\ 3 + 4"e 4A4 + 6£ 4 A 5 + 4e 4 A 6 + e 4 A 7 . 

Here Aq = A" and A£ +1 = 

Consider now the general system (1481) including both the differential equations and the con- 
straints. It follows from Proposition 14.31 that all the gauge symmetries of the differential equations 
alone are also the symmetries of the constraints. Indeed, the distribution £ is tangent to the 
complete constraint surface E, hence (by the Frobenius theorem) the closure £ is also tangent to 
E. This means that not only the vector fields Z a and V respect the constraint surface, but also 
all their iterated commutators do the same: 

(88) XT a = F(X) b a T b VXe£. 

On the other hand, according to Corollary 3 from Proposition I4.3[ the gauge variation S e x l is to 
be spanned by the vectors from the distribution Zy C £. This immediately gives 5 £ T a \^ = 0. 

An important point to stress is that the number of nontrivial gauge transformations can decrease 
when the constraints are taken into account. The matter is that some of the structure functions^ 
R n and U n , which are involved into the gauge transformation ( 149|) . can be trivial ( l23l) with regard 
to the constraints. In other words, some linear combinations of the gauge transformations can 
vanish identically on E. Such transformations are also called trivial. All the gauge symmetry 
transformations are to be considered modulo trivial ones. 

5 As it has been already noticed, all these functions can be chosen to be real analytic. 



26 S.L. LYAKHOVICH AND A. A. SHARAPOV 

To formulate a systematic algorithm for constructing a basis of nontrivial and linearly inde- 
pendent gauge transformations in the presence of constraints we need some algebraic background. 
Given a constrained dynamical system in the complete normal form (I48p . denote by 1Z the ring of 
analytical functions of x l and of a finite number of variables A£. Let X C 1Z denote the principle 
ideal generated by the regular constraints {T a }. For simplicity sake, assume that the ideal X is 
simplqj. Then the quotient 1Z/T is an integrality domain and we can form the field of fractions 
Tx = Ft(TZ/I). The field Tx is a natural substitution for the field of meromorphic functions T 
in the presence of constraints. As a practical matter, it is more convenient to use the following 
equivalent definition of Tx- A meromorphic function / G T is said to be regular if it admits a 
representation / = a/b, where a, b G TZ and b ^ X. The ideal X being simple, all the regular 
meromorphic functions constitute a ring TZj C T. It is easily seen that X is the maximal proper 
ideal of TZj and Tx = TZx/1- Thus Tx is just a subquotient of T and in all practical calculations 
we can replace the elements of the field Tx by their regular representatives in T . 

Now define the n-dimensional vector space Wj over Tx as 



(89) Wx = span^- 



d d 



dx 1 dx r 



Again, we can view Wj as a subquotient of W and represent the vectors of Wj by regular elements 
of the A-distribution W, i.e., those vectors of W whose components are regular meromorphic 
functions of T. 

Since the vector fields V and Z a entering the definition of our dynamical system (1481) are assumed 
to be regular and linearly independent modulo constraints, we can define the m-dimensional 
subspace Zj C Wx as 

(90) Z x = span^ (Z u Z m ) . 

Using the fact that the action of the vector fields V and Z a preserves the constraints T a , one can 
easily see that the M-linear operator D : W — ► W defined by (IMjl induces an M-linear operator in 
the subquotient Wj. We will denote the latter operator by the same symbol D. Similar to the 



An ideal X C 1Z is said to be simple if ab £ X implies either a £ X or b £ X. In the case where 1Z is the ring of 
analytical functions, we have the following criterion of simplicity: If the constraint surface S C W l associated to a 
set of regular constraints T a = is connected, then the principal ideal X — (T a ) is simple. 
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case of unconstrained system we can filter the space Zi by the finite sequence of subspaces 

(91) Z x D D Z?- x) D-OZ^DO, 
where 

(92) = {u G Zj | D p u G Zj U DZj U • ■ ■ U D^Zj} . 

Using this filtration and the Euclidean metric in K n , we can then split Z% in the direct sum of 
subspaces 

(93) Z x = Z\ © Z\ © • • • © Z* , Z v z ~ Z { i ] jZ^ ] . 

Now a straightforward analog of Proposition 14. II for the differential J-j-module Wj states that for 
any w G Z x there exists a sequence of elements Wi,...,w p G Zj such that 

(94) D p w + D p - 1 w l + ■■■ + Dw p „! + w p = 0. 

Let Wo,..., W p G W be regular representatives of the elements w , ...,w p . Then, according to the 
general formulae ( ITDl . the gauge transformations read: 

p—l p—n—l p 

(95) 5 £ X* = E E ^ Dm K-n-m-l , S £ X a = E { Z U V~n , 

n=0 m=0 n=l 

where W 7 ^ = U^Z l a . Let us choose a basis {Z ap } in every subspace To any basis element 
we can associate a gauge transformation of the form ( |95l) and this yields a complete basis of 
undepressible gauge transformations in the presence of constraints. 

We conclude this section by some remarks concerning interpretation of Zy in control theory 
and the theory of gauge systems. Observe that the restriction Zy\-%, being a completely integrable 
distribution of S, endows the constraint surface with the structure of foliation In the 

context of gauge systems, the leaves of this foliation are know as gauge orbits. Two points of the 
constraint surface S are considered equivalent if they belong to the same gauge orbit. Equivalent 
points define the same physical state, so that the space of all physical states of the gauge system 
is identified with the space of leaves E/^-"(E). On the other hand, one of the basic concepts of 
control theory is the notion of attainable set. By definition, a point y G £ belongs to the attainable 
set of a point x G S if one can join y to x by an integral curve of ( 1481) with some fixed functions 
X a (t). In the case where Zy and E are real anaclitic and V G Zy, the so-called orbit theorem [1] 
ensures that the attainable set of a point x G £ coincides with the gauge orbit passing through 
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x. This coincidence is not particularly surprising, since the gauge transformations, involving 
arbitrary functions of time as parameters, allow the (control) functions A a (t) to take on arbitrary 
values at each given instant of time. So, the doctrines of control and gauge theories are in a 
sense complementary: the controllable part of dynamics is non-physical, while the physical part 
is uncontrollable. 

5. Physical observables and an involutive normal form 

By definition, t-local values associated to a dynamical system in the complete form ( 14*8]) are 
functions of the phase-space coordinates x l , undetermined multipliers X a and their derivatives up 
to some finite order. Note that the time derivatives of x l can always be excluded with the help of 
the equations of motion ( )48i) . Therefore, without loss of generality, we can identify the space of t- 
local values with the space of real analytical functions 1Z. Among these functions, there are trivial 
ones (}23l) that vanish on the complete constraint surface S. Two t-local values 0\ and O2 are said 
to be equivalent if their difference is a trivial function. In view of the regularity assumptions, the 
last condition amounts to 

(96) Ox-Oa & 1 -0 2 = F a T a . 

In the previous section, the space of equivalence classes was identified with the quotient 1Z/T, 
where X is the ideal generated by the constraints. 

The infinitesimal gauge transformation (195]) maps any solution of (1481) to another one. Given 
an initial time moment to, the gauge parameters e can be chosen vanishing together with all 
their derivatives involved in the gauge transform. So, the initial data remain the same, while 
the solutions are different. The physical values should evolve in a casual way, i.e., they should 
take the same value on every solution originating from a given initial state. This implies that 
the physical values are to be on-shell invariants of the gauge transformations. As is seen from 
(I95p . the gauge variation of A's starts with the highest time derivative of e's and this derivative 
does not contribute to the gauge transformation of x's. This suggests that the gauge invariant 
t-local values can depend on A's only through the trivial contributions. In other words, each 
t-local physical value can be represented by a function of x l . In view of Proposition 14.31 the gauge 
transformations for the functions of x l are generated by the gauge distribution Zy so that the 
subspace of trivial values is automatically gauge invariant. Thus we are lead to the following 
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definition: A physical observable of a dynamical system brought to the complete normal form f|48|) 
is the equivalence class of a phase-space function 0(x) that is on-shell invariant under the action 
of the gauge distribution Z v , i.e., 

(97) ZO| g = VZeZy. 

For the physical observables, the equations of motion ( j48l) reduce to the form 

(98) = VO, 

where V is the complete drift. The undetermined multipliers X a drop out of these equations as 
the physical observables are invariant under the action of the characteristic distribution Z C Zy. 
Notice that the time derivative of an observable is again an observable because [V, Zy] C Zy and 
the complete constraint surface E is invariant under the action of V and Zy. With initial data 
specified, the unique existence of the solution 0(t), t > t , to equation (1981 follows from two 
facts: (i) the undetermined multipliers are not contained in the equations and (ii) the equivalence 
class of VO is gauge invariant. So, the right hand side of (f9~8j) is the same for any solution x l (t) 
evolving from a given initial state x l (to) = x l Q . All that confirms ones again that the definition of 
the physical observables provides them casual evolution. 

As is seen, the following data are only needed to define the physical observables and their 
time evolution: the phase space U, the complete constraint surface E, the complete drift V, and 
the gauge distribution Zy. The quadruple (U, T^,V,Z V ) can always be unambiguously derived 
from equations of motion in their primary normal form (CQ), ([2]) following the algorithm of the 
previous sections. The converse is not true: different primary normal forms can result in the same 
quadruple (U, E, V, Zy). These differences can be much larger than just the equivalence relations 
( [21]) . ( |22|) for the primary normal form. In particular, the output can be the same, even though the 
algorithm has been applied to dynamical systems with characteristic distributions Z and primary 
constraint surfaces E of different dimensions. With the dynamics in mind of physical observables, 
it seems reasonable to consider two dynamical systems as being equivalent whenever they have 
the same complete constraint surfaces, complete drifts and coinciding gauge distributions^- 

^In control theory, two affine-control systems are called feedback equivalent whenever every solution x(t) of 
equations ([1]) with every given control X(t) of one system coincides with a certain solution of another one, possibly 
with a different control. The feedback equivalence is much more restrictive notion than the definition above: For 
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Among the systems that are physically equivalent to (JT]), (j2J) there is a special one whose 
characteristic distribution, primary constraint set, and primary drift coincide, correspondingly, 
with the gauge distribution, complete constraint set, and the complete drift of the original system: 

(99) x l = V\x) + A a Z£(a;) , T- a {x) = . 

By construction, Z v = span{Z s } is tangent to the complete constraint surface £ = {x G 
U\T a (x) = 0}. This means no compatibility conditions can arise from (1991) . Also, Zy is invo- 
lutive and invariant with respect to the complete drift V. Therefore, the gauge distribution for 
( 1991) coincides with the characteristic distribution. We call equations ( 1991) the involutive normal 
form of local dynamics. The system (1991) describes the observables defined by the same conditions 
( 1971) and having the same evolution law (1981) as the physical observables associated to the original 
equations ([I]), (j2J). As the characteristic distribution is involutive, the gauge transformations for 
( 1991) take quite a simple form (jHJ). 

Let us comment on the involutive normal form (|99l) for a system whose primary normal form 
is the constrained Hamiltonian dynamics PQ, [2]. If equations ([I]), ([2]) follow from the least ac- 
tion principle (H|), then the primary characteristic distribution Z is spanned by the Hamiltonian 
vector fields for the primary constraints and the drift is the Hamiltonian vector field for the 
primary Hamiltonian. As it has been already explained in Section 3, the transverse constraints 
T ± correspond to the second-class constraints, the transverse distribution Z± is spanned by the 
Hamiltonian vector fields for the primary second-class constraints, and Z» is generated by the 
Hamiltonian vector field associated to the primary first-class constraints. Suppose we have no 
transverse constraints. Then the characteristic distribution is tangent to the complete constraint 
surface. From the viewpoint of Dirac's classification, this is the case of a pure first-class system. 
Whenever the Dirac conjecture is truqj, the Hamiltonian vector fields for all the first-class con- 
straints generate the gauge transformations and the involutive form of dynamics (199]) is again 
variational. The corresponding action 

(100) S[x, X) = j (pi^x* - H tot (x, A)) dt , 



feedback equivalent systems all values f(x), both controllable and uncontrollable, have the same time evolutions, 
while for physically equivalent systems, unobservable (i.e., controllable) values may evolve differently. 
8 It is not always true, see [2] for counterexamples. 
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involves the total Hamiltonian H to t — H + X a T a given by the sum of the original Hamiltonian 
([3]) and the linear combination of all the first-class constraints (both primary and secondary) 
with independent Lagrange multipliers. In this case, the gauge distribution is spanned by the 
Hamiltonian vector fields X a = {T a , ■ } for the first-class constraints. 

Whether the Dirac conjecture is true or not, the algorithm of Section 4 shows that the gauge 
distribution is included into span{X s } provided that the primary normal form was variational. 
This does not mean, however, that every Hamiltonian vector field X a should belong to the gauge 
distribution. In case Zy ^ span{X a }, the property conjectured by Dirac does not hold and the 
involutive form of the constrained Hamiltonian dynamics is not variational anymore. Be it as 
it may, the algorithm proposed in this paper will automatically separate the constraints whose 
Hamiltonian vector fields contribute to the gauge transformations from those which Hamiltonian 
vector fields do not. This allows us to systematically identify all the true gauge symmetries for 
any regular, constrained Hamiltonian system, even though the system does not satisfy the Dirac 
conjecture. 

Concluding this section, let us briefly discuss the issue of equipping the involutive dynamics fl99l) 
with a certain Hamiltonian structure when the original equations of motion are not variational. 
The basic idea is that only the algebra of physical observables ( J97l) . not the algebra of all t-local 
values, should be equipped with the Poisson bracket. In other words, it is sufficient to require 
the Jacobi identity to hold only when the Poisson bracket is applied to a pair of observables. 
Besides, the Poisson bracket is to be compatible with the time evolution (|99|) in the sense that 
the time derivative should differentiate the bracket of two observables (not arbitrary functions) by 
the Leibnitz rule. In [5], we introduced a notion of weak Hamiltonian structure, which satisfies 
all the above properties. A similar construction was also studied in [6]. 

Formally, a weak Hamiltonian structure associated to an involutive dynamical system is defined 
by the quadruple (T, Z, V, P), where T = {T a } is a set of constraints, Z = {Z a } are the generators 
of a characteristic(=gauge) distribution, V is a drift, and P = P l WiAdj is a weak Poisson bivector. 
So all the objects are polyvector fields of degree 0, 1, and 2. In terms of the Schouten bracket of 
polyvector fields the defining relations for a weak Hamiltonian structure read 



(101) 



[Z a ,T a } = A b aa T t 



b • 



[Z a , Zp] — BlpZj + T a C\ 



'a 



(102) 



[T a , V] = D b a T b , 



[Z a ,V] = EiZ p + T a F' c 



a , 



<a 
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(103) 



[P,P] 



G a AZ, 



a 



T a H. 



[V, P] = I a A Z a + T a J, 



(104) 



[Ta,P] 



K: AZ a - T b L b a 



[Z a ,P] = M^AZ (3 -T a N' ( 



a 



where A, B,C, N are some polyvector fields. Relations (11 01 1) and H102[) express the fact of 
involutivity of the dynamical system. The first relation in f ll03f) identifies P as a weak Poisson 
bivector, i.e., P satisfies the Jacobi identity modulo constraints and gauge symmetry generators. 
Then the second relation in (11031) tells us that the evolution (1981) generated by V preserves the 
Poisson algebra of physical observables (i.e., V is a weakly Poisson vector field). Relations (11041) 
mean that the Hamiltonian vector fields for the constraints are spanned by the gauge generators 
modulo trivial terms and that the generators are weakly Poisson vector fields. As a result the 
trivial functions constitute the Poisson ideal X C 1Z that makes possible to speak about the Poisson 
algebra of physical observables 1Z/X. 

As it was shown in [5], the weak Hamiltonian structure (I101H1041) admits a nice BRST imbed- 
ding that generalize the usual BFV-BRST formalism for the Hamiltonian systems with first-class 
constraints. Starting with this BRST embedding, it is possible to construct a fully consistent defor- 
mation quantization of a weak Hamiltonian system without any reference to variational principles. 
The output is a weakly associative ^-product inducing an associative quantum multiplication in 
the space of physical observables identified with a certain BRST cohomology. The construction 
essentially relies on the superextension of the formality theorem [6] and may be thought of as a 
generalization of the Kontsevich deformation quantization to the case of non-Hamiltonian gauge 
theories. 

The "odd counterpart" of the weak Hamiltonian structure is known as a Lagrange structure [7J. 
The existence of the latter structure is much less restrictive for the equations of motion than the 
requirement to be variational. Whereas the weak Hamiltonian structure is aimed at the construc- 
tion of ^-product, the Lagrange structure allows one to perform the path-integral quantization of 
a (non-) variational dynamical system. For the variational dynamics, this quantization is shown 
to reduce to the standard BV quantization [TJ, [H]. The analysis of classical gauge symmetries 
performed in Sections 4, 5 may be viewed as a pre-requisite for extending the local BRST co- 
homology techniques [10] from the Lagrangian to non-Lagrangian gauge theories [7J, [8]. Notice, 
however, that to make the technique explicitly covariant in field theoretical context, our analysis 
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should address the issues of locality in multidimensional space. These issues are beyond the scope 
of this work. 

6. Conclusion 

In this paper we propose an algorithm of bringing the general local dynamics to certain normal 
forms, which allow us to identify all gauge symmetries. These normal forms can serve as start- 
ing point for the BRST imbedding and deformation quantization of not necessarily variational 
dynamics. Let us briefly summarize the essentials of the proposed algorithm. 

The algorithm becomes applicable after imposing certain regularity conditions on the local 
equations of motion. Then, we see that the general regular equations can be brought to the 
primary normal form that includes the differential equations (CQ) with undetermined multipliers and 
the phase-space constraints ([2]). In this form, the dynamics are defined by the three ingredients: 
the primary constraints, primary characteristic distribution, and primary drift. In the case of 
variational dynamics, this corresponds to Hamiltonian equations subject to primary constraints. 
The main problem we address in this paper is finding all the gauge symmetry transformations (jSJ) 
for the equations of motion (CQ), (j2J). To this end, we first consider the differential consequences of 
the equations of motion. A basic consistency requirement is the conservation of the constraints in 
time. This results in a multi-step procedure of iterating secondary constraints and determining 
a part of the multipliers. After terminating the procedure, we are left with equations ( 1481) - the 
complete normal form of local dynamics - which assume no further restrictions on the phase-space 
coordinates and/or undetermined multipliers. 

Having the dynamical system brought to the complete normal form, one can go over to finding 
its gauge symmetries. We find that the gauge transformations ( 195|) are generated by the gauge 
distribution Zy. The gauge distribution is the Lie closure of the tangential characteristic distri- 
bution Z supplemented by all its iterated commutators with the drift vector field fj47|) . So, the 
gauge symmetries are explicitly identified in the same terms that define the original system (CQ), 
©• 

Then we turn to the notion of physical observables. These are understood as gauge invariants of 
the dynamics. As the gauge symmetry is generated by the gauge distribution Zy, the observables 
are defined as equivalence classes (1961) of phase-space functions annihilated by Zy. This definition 
( 1971) ensures that the time evolution of observables ( 1981) is casual. It also suggests to consider two 
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dynamical systems as equivalent to each other whenever they have coinciding gauge distributions, 
complete constraint sets and complete drifts. As a result, the physical observables and their time 
evolutions are the same for all the equivalent systems. In every equivalence class of the dynamical 
systems, there is a special representative whose characteristic distribution is involutive, tangent 
to the constraint surface, and is preserved by the drift. In this form, called the involutive normal, 
the gauge transformations take the most simple form (J9j). It is the involutive normal form that 
serves as a starting point for the deformation quantization based on the concept of weak Poisson 
structure fll01til04j> . 
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